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Abstract. I-band CCD images of two large regions of 
the Galactic globular cluster NGC 1261 have been used 
to construct stellar luminosity functions (LF) for 14000 
stars in three annuli from 1'4 from the cluster center to 
the tidal radius. The LFs extend to Mi ~ 8 and tend to 
steepen from the inner to the outer annulus, in agreement 
with the predictions of the multimass King-Michie model 
that we have calculated for this cluster. The LFs have 
been transformed into mass functions. Once corrected for 
mass segregation the global mass function of NGC 1261 
has a slope xq = 0.8 ± 0.5 
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1. Introduction 

One of the basic ingredients of a stellar system is the 
mass distribution of its population of stars. In the case 
of globular clusters (GCs), stellar mass functions (MFs) 
reflect both the initial conditions from the epoch of clus- 
ter formation and the dynamical evolution. In order to 
try to disentangle the effects of evolution from those of 
the initial conditions on the shapes of globular cluster 
MFs we have started a long-term project aimed at get- 
ting deep luminosity functions (LF) in a large sample of 
galactic GCs. The first results of this project have already 
been presented in other papers (c/. Piotto 1993 and Pi- 
otto et al. (1996) for a complete bibliography on recent 
works on GC LFs). Since Scalo (1986) it has been known 
that the MFs of Galactic GCs differ from cluster to clus- 
ter. McClure et al. (1986) proposed that the MF slopes 
depend on the cluster metallicity. By studying a sam- 
ple of 17 Galactic GCs, Capaccioli, Ortolani, & Piotto 
(1991) questioned the hypothesis that the metal content 
could be the main parameter governing the MF shapes; 
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instead, they showed that the position of the parent clus- 
ter within the Galactic gravitational potential is affecting 
the present day MFs (PDMF). By means of a more de- 
tailed statistical analysis on the same sample of Capacci- 
oli et al. (1991), Djorgovski, Piotto, & Capaccioli (1993) 
demonstrated that the MF slopes are determined both 
by the positional variables and by the metallicity. The 
two positional variables are dominant, with similar but 
not equivalent effects: at a given Galactocentric distance 
RgCj clusters with a smaller distance from the galactic 
plane Zgp have shallower MFs, and vice-versa. However, 
they do not determine the MF completely: at a given posi- 
tion, clusters with higher metallicity have shallower MFs. 
This dependence on the positional coordinates has been 
tentatively interpreted as the effect of selective evapora- 
tion of stars, driven by tidal shocks (Capaccioli, Piotto, & 
Stiavelli 1993): heating due to disk shocking acts prefer- 
entially on stars located in the outer regions of the cluster 
and, in combination with mass segregation, it leads to a 
loss of low mass stars, i.e. to a MF flattening. The depen- 
dence on the metallicity could reflect the initial conditions. 
This scenario has been furtherly confirmed by recent HST 
data (Piotto, Cool, & King 1997) which allowed to extend 
the study of GC MFs almost down to the limit of the core 
hydrogen-burning ignition (~ O.IMq). 

In this paper we present I-band stellar luminosity func- 
tions of NGC 1261 =C0310 - 554, a cluster which is 
presently relatively far from the center and the plane of 
the Galaxy. On the basis of the location of NGC 1261 
within the Galaxy, Rgc — 17-8 kpc and Zqp = 12.6 kpc, 
it is very likely that its evolution has been driven more by 
the internal dynamics rather than by strong tidal shocks. 
For this reason, despite its large distance from the Sun, 
NGC 1261 has been chosen (among other clusters) for a 
detailed study of the PDMF, which should be as close as 
the internal dynamical evolution can allow to the initial 
MF (IMF). 

NGC 1261 (azooo = 3''12™15", cJzooo = -55°13'l"; is a 
medium concentration globular cluster (c = 1.27 Trager, 
King and Djorgovski 1995), located at high galactic lati- 
tude, I = 250. 54, b — —52. 13. Its low reddening and small 
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Fig. 1. The two fields observed in 
NGC 1261: the upper one is the 
NTT field, while the lower box 
shows the 2.2ni field 



field star contamination make it a very good target for 
photometric studies. The first CCD color-magnitude dia- 
gram (CMD) for NGC 1261 has been pubhshed by Bolte 
& Marleau (1989). They observed the cluster in B and 
V, down to magnitude V ^ 23. From the apparent mag- 
nitude of the HB and the fitting of the main sequence to 
two subdwarfs, they derived an apparent distance modulus 
(m — M)v ~ 16.05 ± 0.25. By matching their observations 
to the VandenBerg and Bell (1985) isochrones and assum- 
ing E(B-V)=0.02, Bolte and Marleau obtained an age of 
15 Gyr and a metaUicity [M/H]=-1.09. By adopting E(B- 
V)=0.04, the best fitting isochrone gives an age of 15 Gyr 
and a metaUicity [M/H]=-1.27. 

Alcaino et al. (1992) published a BVR photometric 
study of NGC 1261. They found a distance modulus (m — 
M)v = 16.00 ± 0.23, from the apparent magnitude of the 
HB, and, using the same set of isochrones of the previous 
study, an age of 15 ± 2 Gyr, and a reddening of E(B- 
V)= 0.07, assuming a metaUicity [Fe/H]= -1.27. More 
recently, Ferraro et al. (1993) presented a BV CMD and 
a luminosity function for the evolved stars in this cluster. 
From the CMD metaUicity indicators, calibrated with the 
Zinn & West (1984) metaUicity scale, they estimated a 



slightly lower metal content than in the previous studies: 
[Fe/H]= -1.4 ± 0.2. They also detected for the first time 
the 'RGB bump' of the cluster, as a clump of stars in the 
LF of the red giant branch, at magnitude V — 16.70±0.05. 
The paper is organized as follows: in Section 2 we dis- 
cuss our data and the method we followed to obtain the 
photometry. In Section 3 we present the luminosity func- 
tions in three radial annuli. Section 4 concerns the mass 
functions and the correction for the effects of the mass seg- 
regation. Finally, in the last Section the MF of NGC 1261 
is compared with the MFs for other 19 Galactic GCs. 

2. Observation and data reduction 

The database consists of two I-band image sets. The first 
one is a sample of seven 15 niin exposures obtained on 
November 27-30, 1991, with EF0SC2 + CCD # 17 at 
the ESO/MPI 2.2 m telescope of La Silla. These images 
have dimensions of 1024 x 1024 pixels, corresponding to an 
area of 5'8x5'8 on the sky, and cover a region from ~ 1'4 
to ^ 7'0 from the center of the cluster. The second set 
consists of 9 frames, for a total exposure time of 95 min- 
utes, obtained with EMMI + CCD # 34 at the ESO/NTT 
telescope in December 1993. This field cover 9'6 x 8'5 of 



Zoccali et al.: The L.F. and M.F. of NGC 1261 



-1 



-2 



'^^if^^%it^~^ ^^ 



- - r - - ^ 
A 



20 



21 



22 



23 



24 



25 



I. 



r = 8' 





r = 4' r = a'.5 r=r.4 



Fig. 2. This plot shows the difference between the magni- 
tude at which each artificial star was added to the image 
and the magnitude at which it was found, as a function of 
the input magnitude. The graphic refers to the artificial- 
star experiments on the NTT frame. The horizontal and 
vertical dashed lines are the limits of the LF magnitude 
bin (i.e. a star that falls out from the two horizontal lines 
at ±0.250 would have been counted in the previous or in 
the next LF magnitude bin. 

the sky, from ^ 1'2 to ^ 9'5 from the center, as shown 
in Fig. n. All observations have been taken in fairly good 
seeing conditions (< FWHM >'- l^'O). 

The images were bias-corrected, trimmed and flat- 
fielded with routines in MIDAS. The single frames of the 
same field were averaged in order to increase the signal- 
to-noise ratio. Profile-fitting photometry of stars was per- 
formed with DAOPHOT II (Stetson 1987). 

In order to obtain deep luminosity functions it is nec- 
essary to measure the magnitudes of very faint stars with 
good accuracy. This can be dangerous if one has frames 
only in one photometric band, because in this case it is 
very difficult to identify and reject spurious detections. 
For EMMI images we minimized this problem by averag- 
ing separately two sets of frames: a first set of five images 
(say image A) and a second one of four images (image B). 
The stellar photometry on these two images has been done 
independently from each other. The star list used as input 
for the photometry on images A and B was the final list 
of objects resulting from the complete photometry of the 
median image of all the (nine) original frames. This is the 
deepest image, as it has the highest signal-to-noise ratio. 
The final star catalog (used for the LFs) contains only the 
stars measured in both images A and B. The magnitude 



Fig. 3. The position of the stars present in the field ob- 
served with EMMI and the radial annuli in which we 
separated the counts. The inner annulus has boundaries 
1'4 < r < 2'5, the intermediate one has 2'5 < r < 4'0, 
and the outer annulus has 4'0 < r < 8'0 



r = 2',5 



Fig. 4. Same as Fig. |, for the field observed with EF0SC2 



associated to each star is the average of the magnitudes we 
measured in the two images. This method allows to use 
the full information present in our data set and also to 
eliminate most of the spurious detections, as, in general. 
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Table 2. 2.2m Luminosity Functions 
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the noise peaks cannot be present in both images at the 
same position. Indeed, the probabihty that a pixel differs 
by more than 3.5 cr from the average value of the sky is 
0.0004; the probability that the same pixel differs more 
than 3.5 a in both images is 1.6 x 10~^. Since our images 
have dimension of 1700 x 1500 pixels, for a total of 2550000 
pixels, we expect that 0.41 noise peaks should be present 
as spurious detections in our star list. 

The completeness corrections have been estimated by 
standard artificial-star experiments. We performed 10 in- 
dependent experiments per field, adding each time a total 
of 600 stars, with the spatial distribution of a King model 
of the same concentration (c=1.3, Trager et al. 1995) of 
NGC 1261. The finding algorithm adopted to recover and 
measure the artificial stars was the same used for the pho- 
tometry of the original images. As discussed extensively 
by Stetson & Harris (1988) and Drukier et al. (1988), pho- 
tometric errors cause the measured magnitude of an ar- 
tificial star to be different from its input magnitude, just 
as "real" stars may have measured magnitudes different 
from their true magnitudes. 

In general, the measured magnitude of a star is 
brighter than the input magnitude, as shown in Fig. 
(see Stetson & Harris 1988 for a discussion on the origin 
of this phenomenon). A solution to this problem was sug- 
gested by Drukier et al. (1988): using the artificial-star 
data it is possible to set up a two dimensional matrix giv- 
ing, for each input-magnitude bin, the probability that a 
star would appear in each output-magnitude bin. To be 
more explicit, the columns of the matrix represent the in- 
put magnitude bins and the rows represent the output. 
In this way, each matrix element contains the probability 
that a star added in the magnitude bin i (row i) is found 
in the magnitude bin j (column j). This probability is de- 
fined as the ratio of the number of artificial stars found 
in each magnitude bin (rows of the matrix) to the num- 



ber of stars added in a bin (columns of the matrix). The 
inverse of this matrix gives the completeness correction 
factors. The observed LF, once multiplied by this inverse 
matrix, becomes the complete LF, with the two different 
effects of crowding properly corrected: the loss of stars and 
their migration in magnitude (c/. Drukier et al. 1988 for 
the details). In order to evaluate the uncertainty associ- 
ated to each element of the completeness matrix, we per- 
formed nine independent experiments on each field, cre- 
ated nine matrices and inverted each of them separately, 
using the Gauss-Jordan elimination method. By calculat- 
ing the mean value and its standard deviation for each 
element, we were able to determine the errors associated 
to the crowding correction. 

In order to calibrate the instrumental magnitudes to a 
standard system, during each observing night we observed 
27 standard stars at the 2.2m and 18 standards at the 
NTT, in 4 Landolt (1992) standard fields. The calibration 
equations are: 

/,„,t^ ^I^ (0.04 ± 0.01)(y - /) + 2.14 ± 0.01 

for the 2.2m data, and: 

hnstr -I = (0.004 ± 0.002)(y - /) ± 0.60 ± 0.01 

for the NTT data. As we have no V^-band image, we 
adopted a mean color (V — I) — 1.2 for the main- 
sequence stars in the range of magnitudes we were study- 
ing (17< Mi <24). Since the total range in color spanned 
by this region of the main-sequence is A(y — /)=0.8, the 
real color of our stars might be wrong by at most ±0.4 
magnitudes. Using the color term of the two equations 
above, this gives an error of A/ = ±0.02 for the 2.2m data, 
and A/ — ±0.002 for NTT. Considering also the error in 
the zero point, we have a total uncertainty A/ = ±0.03 
for the 2.2m data and A/ = ±0.01 for the NTT data, 
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Fig. 5. Left panel: open circles the radial surface brightness profile of NGC 1261 by Trager et al. (1995). The points 
relative to our counts in the EMMI field (after completeness correction) are represented with filled circles. The fourth 
annulus (which is entirely outside the tidal radius given by Trager et al. ) has been divided into two sub-annuli. Here 
we note that the star counts monotonically decrease till the outer limit of our field, i.e. we never reach the background 
level, even if the last two points refer to stars which are all outside the tidal radius rt = 7'3 given by Trager et al. Right 
panel: comparison between the LF for the annulus outside the published tidal radius (filled circles) corresponding 
to the last two points in the left panel, and the background/foreground counts predicted by Ratnatunga & Bahcall 
(1985) for this region. The observed counts in the right panel are not corrected for completeness. 



negligible in both cases, in view of the fact that the LF 
bin is of 0.5 magnitudes. 

3. Luminosity Functions 

The luminosity function of a system is constructed by sim- 
ply counting the number of stars in successive magnitude 
bins. We separated observed star counts in different radial 
annuli (see Figs. HandH) in order to take into account the 
different effects of crowding at different distances from the 
center of the cluster. Moreover, by determining luminos- 
ity and mass functions at different radial distances we are 
able to study the effects of mass segregation (Section V). 
For both the EMMI and the 2.2 m fields we separated 
the star counts in three annuli: the inner one, from 1'4 
to 2'5, the intermediate one, from 2'5 to 4', and the ex- 
ternal one, from 4' to 8'. The stars in a fourth region of 
the NTT image, with r > 8', have been used to check the 
background/foreground contamination. The dimensions of 
these annuli were chosen in order to have about the same 
statistical sample of stars and to avoid a significant crowd- 
ing gradient inside each one. The latter condition ensured 



that we could apply just one completeness correction for 
each radial luminosity function. 

A peculiarity of NGC 1261 is the fact that it lies on the 
line of sight to a cluster of galaxies. These galaxies were 
identified by DAOPHOT as agglomerates of stars. To solve 
this problem, we visually inspected the images and found 
the coordinates and the sizes of the galaxies directly on 
the frame, and then rejected all the stars within the galaxy 
boxes from the output file of ALLSTAR. 

Another source of contamination is given by the fore- 
ground/background stars that are projected on the clus- 
ter. Normally one can correct this effect by determining 
the luminosity function of a region of sky near the clus- 
ter (in order to have the same population of field stars) 
but far enough to be sure there is no cluster object. Then, 
one has simply to subtract the number of field stars per 
unit area and unit magnitude from the cluster LF. Our 
idea was to consider as a foreground field the region just 
outside the tidal radius of NGC 1261, in the EMMI field. 
As the published value of the tidal radius of NGC 1261 is 
rt = 7'. 3 (Trager et al. 1995), we isolated the stars outside 
8' from the center of the cluster, and determined their lu- 
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minosity function. Unfortunately, a plot of the cumulative 
number of stars per unit area versus the distance from the 
center demonstrates that there is no flattening at the tidal 
radius, i.e. the value of rt seems to be underestimated (or, 
alternatively, there is a halo of cluster stars, beyond the 
tidal radius, as found by Grillmair et al. 1995 for other 
clusters). This fact is confirmed by Fig. in which we 
represent (left panel) the surface brightness of NGC 1261 
as a function of the radial distance. In this plot, the points 
coming from our counts in each annulus have been added 
to the data of Trager et al. (1995). The fourth annulus 
(which is entirely outside the tidal radius given by Trager 
et al. ) has been divided into two sub-annuli. It is clear that 
the surface brightness is still decreasing outside rt , and we 
can not use the LF from the "external" annulus to correct 
for foreground/background contamination. As a further 
evidence, in the right panel of the same figure we compare 
the star counts observed in the region outside 8' with the 
counts predicted by Ratnatunga & Bahcall (1985) model. 
Our counts are overabundant with respect to the expected 
foreground up to a factor of five, well outside the uncer- 
tainty given by Ratnatunga & Bahcall for their model. We 
conclude that all our fields lie entirely inside the boundary 
of the cluster. For this reason, the background/foreground 
correction has been done using the counts by Ratnatunga 
& Bahcall (1985) based on the Galaxy model by Bahcall & 
Soneira (1984). The amount of this correction is, at most, 
10% of the original counts. 

No correction is made here for the possible contamina- 
tion of the LF by fainter, uresolved galaxies. However, a 



simple test was carried out to insure that, would a correc- 
tion be made, the conclusion of this work would have been 
the same. We examined the CMD obtained from two sets 
of deep HST WFPC2 V and I images for the GCs M92 and 
M30, two clusters whose projected position in the Galaxy 
is similar to that of NGC 1261. We measured the number 
of foreground stars and galaxies for these clusters count- 
ing all the objects that lie outside ±5(t from the cluster 
Main Sequence. We then scaled this counts to 1 arcmin^ 
area and subtracted the number of object expected to be 
stars, given by Ratnatunga & Bahcall (1985). The number 
of remaining objects is our best estimate of the number 
of galaxies per arcmin^ that could contaminate our LF. 
Due to their small area, this correction does not affect 
at all the inner and intermediate annuli of our LFs of 
NGC 1261. It would change the last three magnitude bins 
of the outer annuli of both fields. In particular the correc- 
tion would reach up to 5 counts per armin^ in the last bin 
of the NTT frame and 4 counts per arcmin^ in the last 
bin of the 2.2m field. The corresponding MF slopes would 
be xo — 0.6 and xqbs — 1-5, well inside the error bars 
we gave for our adopted values (cfr. section H). The slope 
for the outer MF would make the observed effect of mass 
segregation more similar to what predicted by the model 
(cf. Fig||). However, this correction must be considered 
an upper limit, because the completeness we have in the 
deep HST frames is significantly larger than the complete- 
ness of our ground-based frames for NGC 1261. It results 
quite difficult to estimate the completeness of the back- 
ground galaxy counts. For this reason, and for the very 
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the outer annulus as a consequence of the mass segrega- 
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small statistic (hence large error) we have to deal with, 
when counting foreground objects, we decided not to cor- 
rect our counts. The reader should simply keep in mind 
that our final values are probably an upper limit for the 
outer MF slopes. This slopes are more likely to be com- 
prised between those shown in Fig.QG and the two lower 
limit given above, when applying the foreground galaxy 
correction. 

The radial luminosity functions for each annulus of our 
two fields are listed in Tables 1 and 2. The counts are nor- 
malized per square arcmin. Only counts with a complete- 
ness correction fraction < 2 have been used. Note that, 
where the dominant effect of the crowding is the loss of 
stars, the corrected counts are larger than the raw ones, 
whereas where the dominant effect is the migration of stel- 
lar magnitudes (i.e. in the brighter bins of the LF and/or 
in the outer annuli) the corrected counts are sligtly smaller 
than the original ones. This is due to the combination 
of two effects: the fact that the migration of magnitudes 
moves the stars essentially toward brighter LF bins, and 
the fact that this effect is larger for fainter stars. The errors 
of the incomplete LF represent the standard deviation of 
a Poisson distribution, i.e. the square root of the number 
of stars; the error in the complete counts were obtained 
combining them with the errors in the inverse matrix. The 
outer annulus of the 2.2 m field is the most noisy, because 
it contains a smaller number of stars. Figure shows a 
comparison between the LFs of the corresponding annuli 
of the two fields. There is a good agreement between the 
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Fig. 9. Multimass King-Michie models for a cluster with 
the concentration parameter of NGC 1261. The filled dots 
are the observed values for the slope of the radial MFs. 

two fields, so the LFs used to derive the radial MFs have 
been obtained by adding the LFs of the corresponding 
annuli in the two fields. In the intermediate and outer an- 
nuli we excluded the last magnitude bin of the NTT LF 
as there is not the corresponding bin in the 2.2m data. 
The final LFs for the three annuli are shown in Fig. f\ 
Note how the LF becomes steeper and steeper moving to 
the outer parts of the cluster. As discussed in the next 
Section, this is what we expect as a consequence of mass 
segregation. 

4. Mass Functions 

In order to convert the luminosity function ^{Mj)dMi 
into the mass function ^(rn)dm we need a mass-luminosity 
relation (MLR) appropriate for the Population II stars and 
for the cluster metal abundance. The most recent I-band 
MLR is from Alexander et al. (1997). For NGC 1261 
we used the isochrones corresponding to a metallicity 
[Fc/II]=-1.3 and we performed an interpolation between 
the published isochrones relative to 10 Gyr and 20 Gyr in 
order to obtain an age of 15 Gyr (Alcaino et al. 1992). 
According to the discussion in Section 1, we adopted a 
distance modulus of (m — M) = 16.00 and we neglected 
the absorption which is very small, in view of the small 
reddening and of the fact that we observed in the /-band. 
The mass functions for the three radial annuli are 
shown in Fig. pi. For consistency with the previous work, 
a straight line has been fitted to the data points of both 
fields, in the mass interval 0.5 < m/rriQ < 0.8. The slopes 
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Fig. 10. Trivariate relation by Djorgovski et al. (1993) with the point relative to NGC 1261. The coefRcients have 
been calculated using the new extended set of clusters with available MFs 



of these lines and their errors, obtained by performing a 
weighted least square fit, are written in each panel of the 
figure. 



4-.1. Mass segregation correction 

The slope of the global MP of the stars in NGC 1261 
can be obtained after a correction for mass segregation. 
We used isotropic multiniass King (1966) models following 
Pryor, Smith & McClure (1986) and Pryor et al. (1991). 

Por the multimass models we followed the same recipe 
and used the same code of Pryor et al. (1991). Our only 
change was in the adopted MP, which was assumed to be 
a power law and divided in 7 mass bins from the turn-off 
mass down to 0.1 Mq, plus a mass bin for the remnant 
stars (in total 8 mass bins). The MP slope was varied in the 
range —1.5 < a; < 1.35, finding for each slope the model 
best fitting the cluster radial density profile of Pig. H then 
we calculated the radial variation of x, due to mass segre- 
gation, fitting the local MP in the model, in the same mass 
range of the observed MP: 0.5 < m/rriQ < 0.8. As it can 
be seen in Pig. u, the MP slopes increase from the inner 
to the outer bins, as expected (see also, Pryor et al. 1986). 
Our three values of the MP are consistent with the mass 
segregation predicted by the King models and the global 
MP slope can be estimated to be xq = 0.8 ± 0.5. We also 
determined the total MP, obtained by simply summing 



the number of stars observed in each magnitude interval 
of each annulus. This MP has a slope xqbs = 1.7 ±0.5. It 
is not surprising that this MP is steeper than the global 
one, as it refers to stars located at r > 1'4, i.e. r > 4.4rc, 
and not to the entire stellar population in NGC 1261: this 
is a further evidence of the mass segregation effects. 

The structural parameters of our best fitting dynami- 
cal modelling of NGC 1261 are different from the Trager 
et al values.: Tc = 16" and c — lA while Trager found 
Tc = 27" and c = 1.27. This is mainly because adding our 
star counts to the Trager et al luminosity profile we gave 
a different weight to the outer part of the cluster profile. 
In any case this difference of structural parameters do not 
change our previous results on the global value of x since 
our models are directly best fitted to the data and not 
simply calculated from the structural parameters. 



5. Discussion 

The values of a;o and xqbs have been added to the figures 
with the correlation between the MP slopes of the other 
20 clusters and their position and metallicity as proposed 
by Djorgovski et al. (1993). As shown in Pig. 11^, taking 
into account the uncertainties in the slope, also the MP of 
NGC 1261 follows the correlation proposed by Djorgovski 
et al. (1993). In order to have a quantitative estimate of 
how well the distribution of Pig.|l^ is fitted by a straight 
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line, we calculated the correlation factors, that should be Saviane I., Piotto G., Zaggia S., Fagotto F., Capaccioli M., 

very close to 1 if the distribution is actually a straight line. Aparicio A., 1997, A&A in press 

For the points in the two panels of Fig. p^ we obtained: Scalo J.M., 1986, In "Fundamental of Cosmic Physics", 11, 1 

Stetson P.B. 1987 PASP 99, 101 
r(a;o) = 0.835 and r{xoBs) = 0M3. Stetson P.B. , Harris W.E. 1988 AJ 96, 909 

Trager , King I.R., Djorgovski S. 1995 AJ 109, 218 
We calculated the new coefficients for the trivariate Tyson, J. A. 1988 AJ 96, 1 
formula, including the points for NGC 1261, NGC 3201 VandenBerg D.A., Bell R.A. 1985 ApJS 58, 561 
(Brewer et al. 1993) and NGC 1851 (Saviane et al. 1997) Zinn R., West M. 1984 ApJS 55, 45 

and M55 (Zaggia et al. 1997) using the same method de- Zaggia S.R., Capaccioli M., Piotto C, 1997, A&A in press 
scribed in Djorgovski et al. . The new data suggest that 
the best trivariate relations are: 

Xobs = (2.14 ± 0.85) log Rgc + (0.71 ± 0.36) log Zap - (0.74 ± 0.30) [Fe/H] - 2.66 

for the observed MF slope, and: 
xo = (1.30 ± 0.66) log Rgc + (0.63 ± 0.29) log Zgp - (0.62 ± 0.24) [Fe/H] - 2.03 

for the mass segregation corrected MF slope. 

Acknowledgements. We are grateful to C. Pryor for providing 
us the dynamical code for the mass segregation correction. 

References 

Alcaino C, Liller W., Alvarado F., Wenderoth E., 1992 AJ 

104,1850 
Alexander D.R., Brocato E., Cassisi S., Castellani V., Ciacio 

F., Degl'Innocenti S. 1997 A&A 317,90 
Bahcall J.N., Soneira R.M., 1984 ApJ 55, 67 
Bolte M., Marleau F., 1989, PASP 101, 1088 
Brewer J.D., Fahlman G.G., Richer H.B., Searle L., Thompson 

I., 1993 AJ 105, 2158 
Berghbush P.A., Vandenberg D.A., 1992 ApJ 81, 163 
Capaccioli M., Piotto G., StiaveUi M., 1993 MNRAS 261, 819 
Capaccioli M., Ortolani S., Piotto G., 1991 A&A 244, 298 
Djorgovski S. , Piotto G., Capaccioli M., 1993 AJ 105, 2148 
Druckier G.A., Fahlman G.G., Richer H.B., Vandenberg D.A., 

1988 AJ 95, 1415 
Ferraro F.R., Clementini G., Fusi Pecci F., Vitiello E., Buo- 

nanno R., 1993 MNRAS 264, 273 
King I.R. 1996 AJ 71, 64 
Landolt A., 1992 AJ 104, 340 
McClure R.D., Vandenberg D.A., Smith G.H., Fahlman G.G., 

Richer H.B., Hesser J.E., Harris W.E., Stetson P.B., Bell 

R.A., 1986 ApJL 307,L49 
Pryor C, McClure R.D., Fletcher J.M., Hesser J.E., 1991 AJ 

102, 1026 
Piotto G., 1993, ASP Conf. Ser. 50, 233 
Piotto G., Cool A.M., King I.R., 1997a, lAU Symp. Nr 174, 

"Dynamical Evolution of Star Clusters - Confrontation of 

Theory and Observation", ed. Hut P., Makino J. (Dor- 
drecht: Kluwer Academic Publishers), p. 71 
Piotto G., Cool A.M., King I.R. 1997b AJ 113, 1345 
Pryor C, Smith G.H., Mc Clure R.D. 1986 AJ 92, 1358 
Pryor C, Mc Clure R. D., Fletcher J. M., Hesser J. E. 1991 

AJ 102, 1026 
Ratnatunga U.K., Bahcall J.N., 1985 ApJ 59, 63 
Richer H.B., Fahlman G.G., Buonanno R., Fusi Pecci F., Searle 

L., Thompson I.B. 1991 ApJ 381,147 



